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ABSTRACT. We discuss parabolic versions of Euler's identity 
e !t = cos t + i sin t. 

A purely algebraic approach based on dual numbers is known to produce a very 
trivial relation c £t = 1 + et. Therefore we use a geometric setup of parabolic 
rotations to recover the corresponding non-trivial algebraic framework. Our main 
tool is Mobius transformations which turn out to be closely related to induced 
representations of the group S L 2 ( R ) . 
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I. Introduction: a Parabolic Wheel — an Algebraic Approach 

A mathematical picture of a wheel, which uniformly rotates around its axis, is 
given by the following "model": 

(1.1) x = cost, u=sint, 

where x and y denote the coordinates of a point on the unit distance from the 
axis of rotation. The principal ingredients are sine and cosine functions, which are 
known for more than two thousand years. 

1.1. Complex Numbers. A later invention of complex numbers z = x+ry, i 2 = —1 
allows to write down two identities (1.1) as a single one: 

(1.2) z = cost + isint. 

The next big advance is known as Euler's formula, which expresses trigonometric 
functions through the exponent of an imaginary number: 

(1.3) e [t = cost + isint. 

Thus the geometrical meaning of multiplication by c a is an isometric rotation of 
the plane M 2 , see Fig. 1(E) with the (elliptic) metric given by: 

(1.4) x 2 +y 2 = (x + ru)(x-hj). 
What are possible extensions of this results? 




FIGURE 1 . Rotations of algebraic wheels: elliptic (E), trivial para- 
bolic (Po) and hyperbolic (H). All blue rims of wheels are defined 
by the identity x 2 — i 2 y 2 = 1. Green "spokes" (straight lines from 
the origin to a point on the rims) are "rotated" by multiplication 
by e Lt . 



1.2. Double Numbers. Complex numbers is not the only possible extension of the 
reals. There are other variants of imaginary units, for example 1 e 2 = 1. Replacing 
i by e in (1.3) we get a key to hyperbolic trigonometry: 

(1.5) e et = cosht + e sinht. 

Here expressions x + ey form the algebra of double numbers — the simplest case of 
hypernumbers. Multiplication by e et is a map of double numbers into itself which 
preserves the hyperbolic metric, cf. (1.4): 

(1.6) x 2 -y 2 = (x+ey)(x-ey). 

Geometrically this may be viewed as hyperbolic rotation, see Fig. 1(H), in contrast 
to the elliptic case (1.3). 



Here e is not a real number — a clarification which may be omitted in the case of i 2 = 



-1. 
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1.3. Dual number. To make the picture complete we may wish to add the para- 
bolic case through the imaginary unit 2 of dual numbers defined by e 2 = 0. Since 
e n = for any integer n > 1 we get from Taylor's series of the exponent function 
the following identity: 

(1.7) e £t = 1 + £t. 

Then parabolic rotations associated with c £ 1 acts on dual numbers as follows: 

e £x : a + eh a + z(ax + b). 

This links the parabolic case with the Galilean group [20]. 
Should we conclude, cf . [4, 20], from here that: 

• the parabolic trigonometric functions are trivial: 

(1.8) cospt = l, sinpt = t? 

• the parabolic distance is independent from y, cf. (1.4) and (1.6): 

(1.9) x 2 = (x+ey)(x-ey)7 

• the polar decomposition of a dual number is defined by [ , App. C(30')]: 

v v 

(1.10) u + ev = u(l + e — ), thus |u+£v|=u, arg(u + ev) = — ? 

u u 

• the parabolic wheel looks rectangular, see Fig. l(Po)? 

The analogies (1.3)-(1.5)-(1.7) and (1.4)-(1.6)-(1.9) are quite explicit and widely 
accepted as an ultimate source for parabolic trigonometry [4,15,20]. However we 
will see shortly that there exists a less trivial form as well. 

Remark 1.1. The parabolic imaginary unit £ is a close relative to the infinitesimal 
number £ from non-standard analysis [2, 19]. The former has the property that its 
square is exactly zero, meanwhile the square of the later is almost zero at its own 
scale. 

Remark 1.2. Introduction of double and dual numbers is not as artificial as it may 
looks from the traditional viewpoint, see Rem. 4.1. 

Remark 1.3. In cases when we need to consider simultaneously several imaginary 
units we use i to denote any of i, £, e. 

2. A Parabolic Wheel — a Geometrical Viewpoint 

We make a second attempt to describe parabolic rotations. If multiplication 
(linear transformation) is not sophisticated enough for this we would advance to 
the next level of complexity: linear-fractional. 

2.1. Matrices. Imaginary units do not need to be seen as abstract quantities. We 
may realise them through zero-trace 2x2 matrices as follows: 

with the parabolic £ nicely siting between the elliptic i and hyperbolic e. Then the 
matrix multiplication implies i 2 = —I, £ 2 = ■ I, e 2 = I, where I is the 2x2 identity 
matrix. Correspondingly we have a matrix form of the identities (1.3)-(1.5): 

,_ / t\ / cost sint\ fO t\ /cosht sinhtA 

(2.2) exp(_ t J = (_ sint cost J, exp(^ J^ sinht cosht J- 



'We use different scripts of the epsilon: e denote hyperbolic imaginary unit and £ — parabolic one. 
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However the above pattern is only partially reproduced in the matrix form of (1.7): 

(2.3) exp(° j) = (J \ 

There is also some arbitrariness in our choice of a matrix representation for e, it 
may be equally well given by the lower-triangular form: 

.<-(; jj) implyi n S s).(i ; 

On the first glance this is not a radical difference, however, it does have some 
implications for the Mobius transform from Section 2.3. 

2.2. Cayley Transform. Another matrix form of the identity (1.3) is provided by 
the Cayley transform: 

1 fl —i\ /cost —sin A fl i\ fe lX 

n— it 



(2 5) 

2 \— i 1 J Vsint cost J \ i 1 
where the matrix 

(2 - 6) c -7l(-. 1 

is the Cayley transform from the upper-half plane to the unit disk. It have its 
hyperbolic cousin 

1/1 e N 



(27) C ^V2\-e 1 

which produces a matrix form of (1.5): 

1 / 1 e\ /cosht sinht\ fl — e\ f e et 



^ 2 V~ e 1 J V sinht coshty/ l) ^0 e~ et 

In the parabolic case we use the same pattern as in (2.6) and (2.7): 

1 -£ N 



c £ 



-e 1 



The Cayley transform of matrix (2.3) is: 

1 -e\ /l t\ A e\ _ A + et t \_ fe £t t 
-e 1 Mo 1/ U 1/ ~ V I -eX) ~ \ 



(2.9) 



-et 



This is again not far from the previous identities (2.5) and (2.8), however, the off- 
diagonal (1, 2) -term destroys harmony. 

Remark 2.1. It is not senseless to consider three matrices (2.1), which materialise 
the imaginary units, together. In fact those trace-less matrices form a basis of the 
SI2 Lie algebra of the group SL 2 (M) [ ]. Moreover they are generators of the one- 
parameter subgroups K, N, A correspondingly, which form the Iwasawa decom- 
position SL 2 (R) = ANK of the group SL 2 (R), see Section 4. 

2.3. Mobius maps. The matrix version of Euler's identity from the previous sec- 
tion can be folded back to numbers through the linear-fractional (or Mobius) trans- 
formations. Indeed any 2x2 matrix define a map: 

,r. « n s fa b\ az+ b 
(2.10) , I : z 1 — ► — , where z = u + iv. 



c dj cz+d 

Notably this is a group homomorphism of invertible matrices under multiplication 
(or the group SL 2 (R)) into transformations of conformally completed plane [5,13]. 
More precisely, see [9] or [1 1 ] (an easy-reading), real-valued matrices (2.2) and (2.3) 
act as transformations of the "upper half-plane"; matrices (2.5), (2.8) and (2.9) act 
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FIGURE 2. Rotation of geometric wheels: elliptic (E), two para- 
bolic (P and P') and hyperbolic (H). Blue orbits are level lines for 
the respective moduli. Green straight lines join points with the 
same value of argument and are drawn with the constant "angu- 
lar step" in each case. 



as transformations of the respective "unit disk", see Fig. 2. Those unit disks are 
images of the upper half -planes under respective Cayley transforms [9, § 8]. 

The actions of diagonal matrices from the right-hand side of identities (2.5) 
and (2.8) are straightforward: they are multiplications by c~ 2lt and e~ 2et corre- 
spondingly. Notably the images of the point — i are: 
(2.11) 

3 H \ fe et \ 

it I : — i I— > sin 2t — i cos 2t; I _ e — et J : — e l— * ~ sinh2t — e cosh2t. 

However the parabolic action of matrix (2.9) in formula (2.10) is not such a simple 
one. 

2.3.1. The Upper-Triangular Subgroup. The parabolic version of the relations (2.11) 
with the upper-triangular matrices from N becomes: 

e £t t\ . . ,. a 



(2.12) e _ £t j r-e-t + ett'-l). 

This coincides with the cyclic rotations defined in [20, § 8]. A comparison of this 
result with (2.11) seemingly confirms that ship t = t but suggest a new expression 
for cosp t: 

cosp t = 1 — t 2 , sinp t = t. 
Therefore the parabolic Pythagoras' identity would be: 

(2.13) sinp 2 1 + cosp t = 1, 

which nicely fits in between of the elliptic and hyperbolic versions: 

sin 2 t + cos 2 t = 1, sinh 2 t - cosh 2 t = — 1. 

The identity (2.13) is also less trivial than the version cosp 2 t = 1 from [ ] (see 
also (1.8), (1.9)). 

2.3.2. The Lower-Triangular Subgroup. There is also the second option to define par- 
abolic rotations, it is generated by the lower-triangular variant of the above con- 
struction, cf. (2.4). The important difference now is: the reference point cannot be 
— e since it is a fixed point — as well as any point on the vertical axis. Instead we 
take £ , which is an ideal element (a point at infinity [20, App. C]) since £ is a 
divisor of zero. The proper compactifications by ideal elements for all three cases 
were discussed in [9,13]. 
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We denote the subgroup of lower-triangular matrices by N '. We obtain with it: 



(2.14) 




A comparison with (2.12) shows that this form is obtained by the change t i— > 
t — 1 . The same transformation gives new expressions for parabolic trigonometric 
functions. The parabolic "unit circle" (or cycle [9,20]) is defined by the equation 
x 2 — y — 1 in both cases . However other orbits are different and we will give their 
description in the next Section. 

3. Rebuilding Algebraic Structures from Geometry 

Rotations in elliptic and hyperbolic cases are given by products of complex or 
double numbers correspondingly, however the multiplication of dual numbers 
produces only the trivial parabolic rotation from Fig. l(Po) rather than more in- 
teresting ones from Fig. 2(P) or Fig. 2(P'). Also the coordinate-wise addition of 
vectors on the plane is invariant under elliptic and hyperbolic rotations but is not 
under the parabolic one. Can we find such algebraic operations for vectors which 
will be compatible with parabolic rotations? 

It is common in mathematics to "revert a theorem into a definition" and we will 
use it systematically in this section to recover a compatible algebraic structure. 

3.1. Modulus and Argument. In the elliptic and hyperbolic cases orbits of rota- 
tions are points with the constant norm (modulus): either x 2 + y 2 or x 2 — y 2 . In the 
parabolic case we employ this point of view as well: 

Definition 3.1. Orbits of actions (2.12) and (2.14) are contour line for the following 
functions which we call respective moduli (norms): 

u 2 

(3.1) for N : |u+ ev| =u 2 -v, for N' : |u+ ev| = 7. 

v + 1 

Remark 3.2. (i) The expression | (u, v) | = u 2 — v represents a parabolic distance 
from (0, h) to (u, v), see [9, Lem. 8.4] and is in line with the "parabolic 
Pythagoras' identity" (2.13). 
(ii) Modulus for N ' expresses the parabolic focal length from (0, —1) to (u, v) 
as described in Lem. 8.5]. 

The only straight lines preserved by the both parabolic rotations N and N ' are 
vertical lines, thus we will treat them as "spokes" for parabolic wheels. Elliptic 
spokes in mathematical terms are "points on the complex plane with the same 
argument", thus we again use it for the parabolic definition: 

Definition 3.3. Parabolic arguments are defined as follows: 

(3.2) for N : arg(u+ ev) =u, for N ' : arg'(u+ ev) = — . 

u 

Both Definitions 3.1 and 3.3 possess natural properties with respect to parabolic 
rotations: 

Proposition 3.4. Let w s is a parabolic rotation ofw by angle s = t in (2.12) or s = t _1 
in (2.14). Then: 

|w s | = |w| , arg''' w s = arg''' w + s, 
where primed versions are used for subgroup N '. 



Remark 3.5. Note that in the commonly accepted approach [20, App. C(30')] para- 
bolic modulus and argument are given by expressions (1.10), which are in a sense 
opposite to our conclusions. 
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3.2. Rotation as Multiplication. We revert again theorems into definitions to as- 
sign multiplication. In fact we require an extended version of properties stated in 
Proposition 3.4: 

Definition 3.6. The product of vectors Wi and w 2 is defined by the following two 
conditions: 

(i) arg(wiw 2 ) = argwi + argw 2 ; 

(ii) |wiw 2 | = |wi| ■ |w 2 |. 

We also need a special form of parabolic conjugation. 



Definition 3.7. Parabolic conjugation is given by u + ev = — u + ev. 

Combination of Definitions 3.1, 3.3 and 3.6 uniquely determine expressions for 
products. 

Proposition 3.8. The parabolic product of vectors is defined by formulae: 

(3.3) /orN: (u, v) * (u', v'] = (u + u', (u + u') 2 - (v - u 2 )(v' - u' 2 )); 

(3.4) forN': (u,v) * (u',v') = — ■ — -, — : — ■ — —„ 1 



(u + u') 2 

Although both expressions looks unusual they have many familiar properties: 

Proposition 3.9. Both products (3.3) and (3.4) satisfy to the following conditions: 

(i) They are commutative and associative; 

(ii) The respective rotations (2.12) and (2.14) are given by multiplications with a 
dual number with the unit norm. 

(iii) The product Wiw 2 is invariant under respective rotations (2.12) and (2.14) . 

(iv) The second component of the product ww is jw| 2 . 

3.3. Invariant Linear Algebra. Now we wish to define a linear structure on M 2 
which would be invariant under point multiplication from the previous Subsec- 
tion (and thus under the parabolic rotations, cf. Prop. 3.9(h)). Multiplication by a 
scalar is straightforward (at list for a positive scalar): it should preserve the argu- 
ment and scale the norm of vectors. Thus we have formulae for a > 0: 

(3.5) a-(u,v) = (u, av + u 2 (l - a)) forN, 

(3.6) Q-(u,v) = (u,^-^-l\ forN'. 

On the other hand addition of vectors can be done in several different ways. We 
present two solutions: one is tropical and another — exotic. 

3.3.1. Tropical form. Let us introduce the lexicographic order on K 2 : 

i \ , i m -r j i • f either u < u'; 

(u,v) H (u ,v I if and only if < , , 

1 ' ' y ' ; } [or u = u', v < v'. 

One can define functions min and max of a pairs of points from on R 2 correspond- 
ingly. An addition of two vectors can be defined either as their minimum or max- 
imum. A similar definition is used in tropical mathematics, also known as Maslov 
dequantisation or K min and R max algebras, see [16] for a comprehensive survey. It 
is easy to check that such an addition is distributive with respect to vector mul- 
tiplications (3.5) — (3.6) and consequently is invariant under parabolic rotations. 
Although it looks promising to investigate this framework we do not study it fur- 
ther for now. 
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3.3.2. Exotic form. Addition of vectors for both subgroups N and N ' can be defined 
by the common rules, where subtle differences are hidden within corresponding 
Definitions 3.1 (norms) and 3.3 (arguments). 

Definition 3.10. Parabolic addition of vectors is defined by the following formu- 
lae: 

(,w « arg (/) Wi ■ |wi| (/) + arg (,) w 2 ■ \W2\_ 



(3.7) arg l ' J (wi+w 2 



|wi + w 2 | (/) 



(3.8) \w 1 +w 2 \ U] = | Wl | ( " + | W2 | ( ", 

primed versions are used for the subgroup N '. 

The rule for the norm of sum (3.8) may looks too trivial at a first glance. We 
should say in its defence that it nicely sits in between of the elliptic jw + w'| 
|w| + |w'| and hyperbolic |w + w'\ ^ |w| + |w'| inequalities for norms. 

Both formulae (3.7)-(3.8) together uniquely define explicit expressions for addi- 
tions of vectors. Although those expressions are rather cumbersome and not really 
much needed. Instead we list properties of this operations: 

Proposition 3.11. Vector additions for subgroups N and N ' defined by (3.7)-(3.8) satisfy 
to the following conditions: 

(i) They are commutative and associative. 

(ii) They are distributive for multiplications (3.3) and (3.4); consequently: 

(iii) They are parabolic rotationally invariant; 

(iv) They are distributive in both ways for the scalar multiplications (3.5) and (3.6) 
respectively: 

a ■ (wi + w 2 ) = a ■ Wi + a ■ w 2 , (a + b) • w = a ■ w + b ■ w. 

To complete the construction we need to define the zero vector and inverse. 

Proposition 3.12. (N) The zero vector is (0,0) and consequently the inverse of (u,v) 
is (u, 2u 2 — v). 

(N ') The zero vector is (00, —1) and consequently the inverse of (u, v) is (u, — v — 2). 

Consequently we can check that scalar multiplications by negative reals are 
given by the same identities (3.5) and (3.6) as for positive ones. 

3.3.3. The Real and Imaginary Parts. Having the vector addition at our hands we 
may wish to define the real and imaginary parts compatible with it. We can start 
from the familiar formulae ^(w + w) and i(w — w). While such a real part has 
a reasonable value (0, |w|) (the subgroup N case), the imaginary part suffers from 
the malformed denominator in (3.7) due to the fact |w| = |w|. 

This should not be seen as a defect of the exotic addition. A moment of reflec- 
tion reveals that "purely real" dual numbers are naturally defined as having zero 
argument, e.g. the vertical axis for the subgroup N . This seems to agree well with 
the above real part. 

However what are "purely imaginary" dual numbers? The horizontal axis can 
be the first suggestion, but all its points have different arguments. One may still 
prefer to choose that the number (1,0) would be purely imaginary (since (0, 1) 
is purely real) and all imaginary numbers have the same argument. So purely 
imaginary numbers 

Definition 3.13. For both subgroups N and N ': 

(i) purely real numbers are defined by the condition argw = 0. 

(ii) purely imaginary numbers are defined by the condition argw ~ 1. 
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Since arguments of 9hv and 3w are fixed by this Definition we need only to find 
their moduli. 

Proposition 3.14. The natural condition w = D^w + 3w together with Defn. 3.13 
uniquely define $Rw and Jw. It is also determined by the identities: 

lEHwl = (1 — argw) |w| , |Jw| = argw |w| . 

The explicit formulae for the real and imaginary parts in the cases of both sub- 
group N and N ' can be found in App. A. 

3.3.4. Linearisation of the exotic form. Some useful information can be obtained from 
the transformation between the parabolic unit disk and its linearised model. In 
such linearised coordinates (a, b) the addition (3.7)-(3.8) is done in the usual coor- 
dinate-wise manner: (a, b) + (a', b') = (a + a', b + b'). 

To this end we calculate the value of (u, v) = a-(ui, Vi)+b-(u2,v 2 ) for (ui,Vi) 
1 . 0) and (u2, V2) = (—1, 0). For the subgroups N the transform is given by: 

(3.9) 

a-b (a-b) 2 . u 2 -v_ . , u 2 -v_ , 

u=— r, v = -a + b), a=— — 1 + u, b = — — 1-u. 

a + b (a + o) z 2 2 

For the subgroup N ' such a transformation is: 

, aim a + b (a + b) u(u+l) , u(u-l) 

(3.10) u = -, v = r-Tj — 1] a 



a-b' (a-b) 2 ' 2(v + l)' 2(v+l) 

We also note that both norms (3.1) have exactly the same value a + b in the respec- 
tive (a, b) -coordinates. 

Remark 3.15. The irrelevance of the standard linear structure for parabolic rotations 
manifests itself in many different ways, e.g. in an apparent "non-conformality" of 
lengths from parabolic foci, e.g. with the parameter cr = in [9, Prop. 5.10.(iii)]. 
An adjustment of notions to the proper framework restores the clear picture. 

The initial definition of conformality [9, Defn. 5.9] considered the usual limit 
y ' — > y along a straight line, i.e. "spoke" in terms of Fig. 1. This is justified in the 
elliptic and hyperbolic cases. However in the parabolic setting the right "spokes" 
are vertical lines, see Fig. 2, so the limit should be taken along them [ , Prop. 5.11]. 

4. Induced Representations as a Source of Imaginary Units 

As we already mentioned in Rem. 2.1 all three matrix exponents (2.2)-(2.3) are 
one-parameter subgroups of the group SL 2 (R) — the group of 2 x 2 matrices with 
unit determinant. Moreover any one-parameter subgroup of SL 2 (R) is a conjugate 
to either of subgroup A, N or K, see Rem. 2.1 for their descriptions. Thus our 
consideration may be applied for construction of induced representations of SL 2 (R) 
in an extended meaning. 

The general scheme of induced representations is as follows, see [6, § 13.2], 
[8, § 3.1]. We denote SL 2 (M) by G and let H be its subgroup. Let CI = G/H be the 
corresponding homogeneous space and s : CI — > G be a continuous function [6, 
§ 13.2] which is a left inverse to the natural projection G — > G/H. In our case we 
choose: 

1 (v u 



(4.1) S: ( U)V )^_I i j j ( U! v)er,v>0. 

Then any g e G has a unique decomposition of the form g = s(u))h where w E CI 
and h G H. We will write: 

(4.2) cu = s- 1 (g), h = r(g):=(s- 1 (g))" 1 g. 
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Note that CI is a left homogeneous space with the G-action defined in terms of s as 
follows: 

(4.3) g : w h-> g • w = s _1 (g _1 * s(w)], 

where * is the multiplication on G . 

Let x : H — > M 2 be a "unitary character" of H in some generalised sense illus- 
trated bellow. Then it induces a "unitary" representation of G, which is very close 
to induced representations in the sense of Mackey [6, § 13.2]. This representation 
has the canonical realisation p in the space L 2 (0) of square integrable K 2 -valued 
functions. It is given by the formula (cf. [6, § 13.2.(7)-(9)]): 

(4.4) [p x (g)f](a))=xo(r(g- 1 *s(a ) )))f(g-c U ), xo(H) =x(H) ( ^j^p ) ' » 

where g e G, m e Q, li 6 H and r:G^H, s:£l^G are maps defined above; 
* denotes multiplication on G and • denotes the action (4.3) of G on D from the 
left. 

4.1. Induction from K. This is the most traditional case in the representation the- 
ory. The action (4.3) takes the form: 

'a b\ / (au + b)(cu + d) + cav 2 v 



c dj-' ' ' \ (cu+d) 2 + (cv) 2 ' (cu+d) 2 + (cv) 2 . 

Obviously it preserves the upper-half plane v > 0. Moreover with the help of the 
imaginary unit i 2 = — 1 it can be naturally represented as a Mobius transformation: 

a b\ aw + b 

. :wh , w = u + iv. 

c d) cw + d 

Thus it is justified in this particular case to look for complex-valued characters of 
K. They are parametrised by an integer n € Z: 

/ cos t — sin 1 1 



p-n^ht) = e mt = (cost + isint) n , where h t = 
We can also calculate that: 

r(g _1 * s(w)) = 



\ sin t cos t 



/cu + d — cv 



^(cu + d) 2 + (cv) 2 V cv cu + d 



II - — 

Taking into account the identity ^ = (-^) 2 we obtain such a realisation of (4.4): 

[Pn(fl)fl(w) = (c^) 2 f (^Tl)' Where ^=(c d)' W=U + iv - 
4.2. Induction from A. In this case the action (4.3) takes the form: 
'a b\ / (au + b)(cu + d) — cav 2 v 



k c dj-' ' ' \ (cu+d) 2 -(cv) 2 ' (cu + d) 2 -(cv) 

This time the map does not preserve the upper-half plane v > 0: the sign of (cu + 
d) 2 — (cv) 2 is not determined. To express this map as a Mobius transformation we 
require the double numbers imaginary unit e 2 = 1: 



a b\ aw + b 

, : w i— ► , w = u + ev. 

c d cw+d 



Remark 4.1. As we can see now the double numbers naturally appear in relation 
with the group SL 2 (R) and thus their introduction in § 1.2 was not "a purely gen- 
eralistic attempt", cf. [17, p. 4]. The same is true for dual numbers as can be seen 
in the next subsection. 
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Under such conditions it does not make much sense to look for a complex- 
valued characters of A. Instead we will take double number valued characters 
which are parametrised by a real number cr: 

p ff (h t ) = e e<rt = (cosht + e sinhtr, where h t = ^ 

We can also calculate that 

, _i , v, 1 f cu+ d cv 

r(g *s(u)J) 



^(cu + d) 2 - (cv) 2 V cv cu + d 
Thus the formula (4.4) becomes: 

[p g (g)f](w) = ( CW + ^ ) f ( aW + \ ) , where g- 1 ^!^ y\ ,w = u+ev. 



cw + dy \ cw + d y ' \c d 

4.3. Induction from N. We consider here the lower-triangular matrices forming 
the subgroup N '. The action (4.3) takes now the form: 

a b\ , , / au + b 
u,v) 



dj ' ' \cu + d' (cu+d) 2 

This map preserves the upper-half plane v > as the elliptic case of K. To express 
this map as a Mobius transformation we require the dual numbers imaginary unit 

e 2 = 0: 

'a b\ aw ■ 



,:wh , w = u + ev. 

c d) cw + d 

Similarly to the previous hyperbolic case we would look not for complex-valued 

characters but rather use parabolic rotations described in Section 2.3. A distinction 

from the hyperbolic case is that they are not given by multiplication of double 

numbers. Such a "character" parametrised by a real number k is defined by: 

p K (ht) : w ~ (1 + £Kt)W + Kt = (l+2£Kt)w+(Kt+£K 2 t 2 ), where h t = (} J 

1 — £Kt Vt 1 



Furthermore we calculate that 



r(g 1 * s(tu)] = _ 

v cu+d 

Thus the formula (4.4) has the following realisation: 



1 

VC 1 



2kvc \ /aw + b\ kvc (kvc) 2 
[p K (g)f](w) = 1 - £ — - f 



cu+d/ \ cw + dy cu + d (cu + d) 2 

where g _1 ~ d) ' w ~ u + £v - 

The vector space of functions, where this representation acts, should be also con- 
sidered with linear operations defined in § 3.3.2. 

These three examples will be used to build the corresponding versions of the 
Cauchy integral formula along the Erlangen Program at Large outlined in [8,9]. 

Acknowledgments 

I am grateful to Prof. S.L. Blyumin for pointing out the paper [ ) 6] to my atten- 
tion. Anastasia Kisil read the manuscript and made useful suggestions. 



12 



VLADIMIR V.KISIL 



References 

[I] Christian Bauer, Alexander Frink, Richard Kreckel, and Jens Vollinga, GiNaC is Not a CAS. 

[2] Martin Davis, Applied nonstandard analysis, Wiley-Interscience [John Wiley & Sons], New York, 

1977. Pure and Applied Mathematics.MR # MR0505473 (58 #21590) 
[3] GNU, General Public License (GPL), version 3, Free Software Foundation, Inc., 59 Temple Place - 

Suite 330, Boston, MA 02111-1307, USA, 29 June 2007. http : //www . gnu . org/licenses/gpl . html. 
[4] Francisco J. Herranz, Ramon Ortega, and Mariano Santander, Trigonometry of spacetimes: a new self- 
dual approach to a curvature/signature {independent trigonometry, J. Phys. A 33 (2000), no. 24, 4525- 

4551. E-print: arXiv:math-ph/9910041.MR #MR1768742 (2001k:53099) 
[5] Francisco J. Herranz and Mariano Santander, Conformal compactification of spacetimes, J. Phys. A 35 

(2002), no. 31, 6619-6629. E-print: arXiv:math-ph/0110019.MR # MR1928852 (2004b:53123) 
[6] A. A. Kirillov, Elements of the theory of representations, Springer- Verlag, Berlin, 1976. Translated 

from the Russian by Edwin Hewitt, Grundlehren der Mathematischen Wissenschaften, Band 

220.MR # 54 #447 

[7] Vladimir V. Kisil, Erlangen program at large — 2: Inventing a wheel. The parabolic one. Preprint LEEDS- 

MATH-PURE-2007-07. E-print: arXiv: 0707. 4024. 
[8] , Analysis in R 1 ' 1 or the principal function theory, Complex Variables Theory Appl. 40 (1999), 

no. 2,93-118. E-print: arXiv:funct-an/9712003. MR # 2000k:30078. 
[9] , Erlangen program at large — 1: Geometry of invariants (2005). E-print: 

arXiv : math . CV/0512416. preprint LEEDS-MATH-PURE-2005-28. 
[10] , An example of Clifford algebras calculations with GiNaC, Adv. in Appl. Clifford Algebras 15 

(2005), no. 2, 239-269. E-print: arXiv : cs . MS/0410044. MR # 2007a:15042. 

[II] , Erlangen program at large — 0: Starting with SL(2, R) group (2006). E-print: 

arXi v : math . GM/0607387. preprint LEEDS-MATH-PURE-2006-11. 

[12] , Pillmore-Springer-Cnops constructions implemented in GiNaC, Adv. Appl. Clifford Algebr. 

17 (2007), no. 1, 59-70. Updated full text and source files: E-print: arXiv : cs . MS/0512073. 
[13] , Two-dimensional conformal models of space-time and their compactification, J. Math. Phys 48 

(2007), no. 073506. E-print: arXiv:math-ph/0611053. 
[14] Serge Lang, SL-2 (R), Graduate Texts in Mathematics, vol. 105, Springer- Verlag, New York, 1985. 

Reprint of the 1975 edition.MR # 803508 (86j:22018) 
[15] Lavrent'ev, M. A. and Shabat, B. V, Problemy gidrodinamiki i ikh matematicheskie modeli. (russian) 

[problems of hydrodynamics and their mathematical models], Second, Izdat. "Nauka", Moscow, 1977. 

MR # 56:17392.MR # 56 #17392 
[16] G. L. Litvinov, The Maslov dequantization, and idempotent and tropical mathematics: a brief introduction, 

Zap. Nauchn. Sem. S.-Peterburg. Otdel. Mat. Inst. Steklov. (POMI) 326 (2005), no. Teor. Predst. Din. 

Sist. Komb. i Algoritm. Metody. 13, 145-182, 282. E-print: arXiv:math/0507014.MR # MR2183219 

(2006i:46104) 

[17] L. S. Pontryagin, Obobshcheniya chisel, BH6jinoTeiiKa "KBaH-r" [Library "Kvant"], vol. 54, 
"Nauka", Moscow, 1986.MR # MR886479 (88c:00005) 

[18] Norman Ramsey, Noweb — a simple, extensible tool for literate programming. 
http: //www. eecs .harvard. edu/"nr /noweb/. 

[19] V. A. Uspenskii, Chto takoe nestandartnyi analizl , "Nauka", Moscow, 1987. With an appendix by V. 
G. Kanovei.MR # MR913941 (88m:26028) 

[20] I. M. Yaglom, A simple non-Euclidean geometry and its physical basis, Springer- Verlag, New York, 
1979. An elementary account of Galilean geometry and the Galilean principle of relativity, Heidel- 
berg Science Library, Translated from the Russian by Abe Shenitzer, With the editorial assistance 
of Basil Gordon.MR # MR520230 (80c:51007) 

Appendix A. Output of Symbolic Calculations 

Here are the results of our symbolic calculations. The source code can be ob- 
tained from this paper [7] source at http://arXiv.org. It uses Clifford algebra 
facilities [ ] of the GiNaC library [1]. The source code is written in noweb [18] 
literature programming environment. 

Cayley of the matrix x: 

Rotation by x: ( u + x x z + 2ux + v ) 
Rotation of (un, v ) by x: ( x — 1 + x 2 ) 
Parabolic norm: — v + u 2 
Real number x as a dual number: ( — x 



1 e°x 

-fe 1 *e°)x + l 
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Product: ( u + u' -u 2 u' 2 + u' 2 + u' 2 v - v'v + 2uu' + uV + u 2 ) 
Product by a scalar: ( u qv + u 2 — u 2 a ) 
Real part: ( -uv + v — u 2 + u 3 ) 
Imag part: ( 1 1 + uv — u 3 ) 

Lin comb of two vectors a»(l, 0)+b»(-l, 0): ( ^ - 3ab2+2ab + b ( 3 - + b b 2 - a2 + 3a2b - 

P is the sum Re(P) and Im(P): true 

The real part of a real dual number is itself: true 

norm is invariant: true 

Product is invariant: true 

Product is norm squared: true 

Product (u,v) * (uo,v ) is (u,v): true 

Add is commutative: true 

Add is associative: true 

S-mult commutative: true 

S-mult associative: true 

S-mult distributive 1: true 

S-mult distributive 2: true 

Product is symmetric (commutative): true 

Prod is associative: true 

Product is distributive: true 



Cayley of the matrix x: 
Rotation by x: 



(e 1 * e°)x + 1 



u 2 x 2 — 2ux- v 



-1+ux l+u 2 x 2 — 2ux 
1 -1+x 2 



Rotation of (uo, vo) by x: 
Parabolic norm: tt- 

l+v 

Real number x as a dual number: ( oo 



Product: ^ 



— l+u' 2 — v' — v'v+2uu'— v+ix 2 



u+u' (u+u'r 

— 1+a— v 



Product by a scalar: ( u 
Realpart:(oo ^gg^ 
Imag part: ( 1 ^ 



u 



2ab-b 2 -a 2 + b 
(a-b) 2 



Lin comb of two vectors a*(l, 0)+b*(-l, 0): { 

P is the sum Re(P) and Im(P): true 
The real part of a real dual number is itself: true 
norm is invariant: true 
Product is invariant: true 
Product is norm squared: true 
Product (u,v) * (uo,v ) is (u, v): true 
Add is commutative: true 
Add is associative: true 
S-mult commutative: true 
S-mult associative: true 
S-mult distributive 1: true 
S-mult distributive 2: true 
Product is symmetric (commutative): true 
Prod is associative: true 
Product is distributive: true 
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Elliptic case of induced representations 

d 2 dc 

map r(M) 



map s 
character: 




cv(d+uc) 



2udc + c 2 v 2 + d 2 +u 2 c 2 2udc + c 2 v 2 + d 2 +u 2 c 2 I 

cvjd+uc) (d+uc) 2 I 

2udc + c 2 v 2 + d 2 +u 2 c 2 2udc + c 2 v 2 + d 2 +u 2 c 2 / 

ucb+db+uad + acv 2 +u 2 ac adv— cvb 

2udc + c 2 v 2 + d 2 +u 2 c 2 2udc + c 2 v 2 + d 2 4 



Moebius map: 
Moebius map is given by the imaginary unit: true 



Parabolic (N') case of induced representations 

mapr(M): ^ \ ° 

( - b 
maps-^M): I jj d 



character: 




1 



Moebius map: ^ -3+ 



cv 
d+uc 

b adv— cvb 



Moebius map is given by the imaginary unit: true 
Hyperbolic case of induced representations 




d 2 


dc \ 


d 2 -c 2 
dc 


d 2 -c 2 
d 2 


d 2 -c 2 


d 2 -c 2 / 




c 2 b + c-d 2 l 


a 


d 2 
d 2 -c 2 
d 


2udc- 


l-d 2 +u 2 c 2 



cv(d+uc) 



2udc — c 2 v 2 + d 2 +u 2 c 2 2udc — c 2 v 2 + d 2 +u 2 c 2 
cv(d+uc) ( d+uc) 



2udc-c 2 v 2 + d 2 +u 2 c 2 2udc — c 2 v 2 + d 2 +u 2 c 2 

Moebius map: ( ^t^ a & ^X> C 2udc^ 
Moebius map is given by the imaginary unit: true 



Appendix B. Program for Symbolic Calculations 

This is a documentation for our symbolic calculations. You can obtain the pro- 
gram itself from the source files of this paper [7] at arXiv . org; fflgX compilation 
of it will produces the file parab-rotation.nw in the current directory. This is a 
noweb [18] code of the program. It uses Clifford algebra facilities [12] of the GiNaC 
library [1]. 

This piece of software is licensed under GNU General Public License (Version 
3, 29 June 2007) [3]. 

B.l. Class dual_number. 

B.l.l. Public Methods. 

A dual number can be created simply by listing its two components. 

14 (Public methods 14)= (28c) 15a > 

dual_number(const ex & a, const ex & b); 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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Alternatively you can provide a 1 x 2 or 2 x 1 matrix, a list, another dual number 
P or a complex expression with a non-zero imaginary part to give two compo- 
nents. If P does not have two components and is a real-valued expression, it will 
be embedded into dual numbers with zero argument and norm equal to P. 

15a (Public methods 14)+= (28c) < 14 15b > 

dual_number(const ex & P); 

Uses dual-number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

We can also obtain the module and argument of a dual_number. 

15b (Public methods 14)+= (28c) <15a 15c> 

ex argQ const; 
ex normQ const; 



We define the conjugate of a dual_number by u + ev = — u + ev. 

15c (Public methods 14)+= (28c) <15b 15d> 

ex conjugate^ const { return dual_number(-u_comp, vjcomp); } 



Uses dual-number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 



Real part may be defined out of the formula SH(u,v) = |((u,v) + (u, v)), but it 
simply reduces to the value of norm for the dual_number. 

15d (Public methods 14)+= (28c) «15c 15e> 

ex reaLpartQ const; 
ex imagjpartQ const; 

Negative of a dual_number and its power. 

15e (Public methods 14)+= (28c) <15d 15f> 

dual_number negQ const 

{ return dual_number(M_comp, -vxomp + (issubgroupjsl? 2*poiv(uxomp, 2) : -2)); > 
dual_number poit>er(const ex & e) const; 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

We can also convert a dual_number to a matrix. 

15f (Public methods 14)+= (28c) <15e 15g> 

matrix tojnatrixQ const { return matrix (1, 2, \st{uxomp, vxomp)); } 

We define the rule for parabolic norm of a sum, see (3.8). 

15g (Public methods 14)+= (28c) <15f 15h> 

ex fldcLnorffxs(const dual_number & P) const { return (normQ+P. normQ). normalQ; } 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

Algebraic operations are defined for dual_numbers in a way described in § 3.3.2. 
The standard C++ operators +, -, *, -r- will be overloaded later in order to permit 
natural expressions with dual_numbers. 

15h (Public methods 14)+= (28c) <15g 

dual_number add(const duaLnumber & a) const; 

dual_number sufo(const dual_number & a) const { return add(a.neg()); } 

dual_number mw/(const dual_number & a) const; 

dual_number mw/(const ex & a) const { return mi(/(dual_number(fl)); > 



Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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B.2. Algebraic Subroutines. We need a couple of global variables which help to 
write uniformly algebraic rules for both cases of subgroups N and N'. 
Firstly, we need to consider separately cases of subgroup N and N', the following 
global variable keeps track on it. 
16a (N-Nprime separation 16a) = (24b) 16b > 

bool issubgroupjsl; 

In the case of the subgroup N' the reference point lies at infinity the following 
realsymbol variable represents it in the symbolic calculations. 

16b (N-Nprime separation 16a) += (24b) < 16a 16c > 

realsymbol In/("Inf ", "Winfty"); 

We define the "zero angle": for the subgroup N it is 0, for N' — oo. 
16c (N-Nprime separation 16a) += (24b) < 16b 

#define ArgO (is_subgroup_N ? ex(0) : ex(Inf)) 

Defines: 

ArgO, used in chunks 21b, 25f, and 30b. 
Here is the set of algebraic procedures representing definitions made in this paper. 

B.2.1. Argument and u. In the case of N the value of u is simply the argument. 
16d (Algebraic procedures 16d) = (24c) 16e> 

ex ufromjirgiconst ex & a) { 
if (issubgroiipJST) 
return a; 

In the case of N ' the value of u is the inverse to the argument, and we need to treat 
properly the case of zero. . . 
16e (Algebraic procedures 16d)+= (24c) <16d 16f> 

else { 

if (a.normalQ.is-zeroQ) 
return Inf; 

. . . and oo. We try to replace — by 0. 

16f (Algebraic procedures 16d)+= (24c) <16e 17e> 

else 

try-C 

realsymbol f("t"); 

return pow(a.subs(Inf= pow(t, -1)), -T).normal().subs(t = 0).normal(); 
} catch (stdv.exception &p) { 
return pow(a, -1); 

} 

} 

} 



B.2.2. Argument of a point. The opposite task (finding argument of a point) is solved 
similarly. 

16g (Dual number class further implementation 16g)= (24c) 17a > 

ex dual_number: :arg() const { 
if (issubgroupjsl) 
return tijcomp; 



Uses dual-number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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Again in the case of N' we need to consider cases of 0. . . 
17a (Dual number class further implementation 16g}+= (24c) < 16g 17b > 

else { 

if (u_comp.normal().is-zero()) 
return Inf; 

. . . and oo. We try to replace by 0. 

17b (Dual number class further implementation 16g}+= (24c) < 17a 17c > 

else 

try-C 

realsymbol f("t"); 

return poiv(iijcomp.subs(Inf= pow(t, -1)), -l).normal().subs(t = 0).normalQ; 
} catch (st d:\exception &p) { 
return pow(uj:omp, -1); 

} 



B.2.3. Norm. The corresponding value of the parabolic norm is calculated by the 
formulae (3.1). 

17c (Dual number class further implementation 16g)+= (24c) <17b 17d> 

ex dual_number::nor»/() const { 
if (is subgroup _N) 

return pow(ujcomp, 2)-vjcomp; 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

The case of subgroup N ' require treatment of infinity. 

17d (Dual number class further implementation 16g)+= (24c) < 17c 18b > 

else { 

if (ujcomp.isjzero()) { 

if ((v_comp+l).is_zeroQ) 
return 1; 

else 

return pow{lnf, 2)~(vjcomp+l); 

} else 

return {pow(uxomp , 2)~{vjcomp+\)).normal(); 

> 



B.2.4. The value of v from the argument and norm. We oftenly need to find values of 
v such that for a given value of argument A point (A, v) will have a given norm. 

17e (Algebraic procedures 16d)+= (24c) <16f 18a > 

ex v-fromjiorm(coT\st ex & u, const ex & n) { 
realsymbol 1{"V); 
if (issubgroupjsl) 

return Zsofoe(dual_number(w, l).norm() = n, l).normal(); 

else 

return /sofoe(poiy(dual_number(w, l).norm(), -1) = pow(n, -1), l).normal(); 

y 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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18a (Algebraic procedures 16d)+= (24c) <17e 18c > 

duaLnumber zerojluaLnnmberQ { 

return (issubgroupjsl ? dual_number(0, 0) : dual_number(7w/, -1)); 

} 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

18b (Dual number class further implementation 16g)+= (24c) <17d 18d> 

bool duaLnumber::zs_zero() const { 
return isxqual{zerojdualjiumberQ); 

} 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

B.2.5. Real and Imaginary Parts. See § 3.3.3 for a discussion of the real and imagi- 
nary parts of dual numbers. 

18c (Algebraic procedures 16d)+= (24c) < 18a 19a > 

ex dual_number: ;real-part() const { 

return dnfromjirgjnod(0, (l-arg())*norm()); 

} 

ex dual_number::zmflg_part() const { 

return dnfromjirgjnod(l, arg()*normQ); 

} 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

B.2.6. Product of Two Points. We define now the product of two points according to 
the Definition 3.6. We also include a multiplication by a scalar: if a factor is a scalar 
it is replaced by a vector with the zero argument and norm equal to the scalar. 

18d (Dual number class further implementation 16g)+= (24c) <18b 18e> 

dual_number dual_number::mu/(const dual_number & P) const { 

ex u=nfromjirg(arg()+P.arg()).normal(); 

return dual_number(i;, V-fromjnorm(u, norm()*P.norm())); 

} 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

B.2.7. Vector Addition of Two Points. The sum is calculated from the expression (3.7). 

18e (Dual number class further implementation 16g)+= (24c) <18d 19b > 

duaLnumber duaLnumber: :«dd(const duaLnumber & a) const { 

ex norms = addjnorms(a); 
if (norms.normalQ.isjzeroQ) 

return zerojiuaLnumberQ; 
else { 

ex us=ufromjirg((argQ*norm()+a.arg()*a.norm()) -f- norms). normalQ; 
return dual_number(us, v_fromjiorm{us, norms)); 

y 

} 



Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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19a (Algebraic procedures 16d)+= (24c) <18c 

duaLnumber dn_fromjirgjnod(const ex & a, const ex & n) { 
ex us=u_fromjirg{a).normal§; 
return dual_number(ws, v_fromjwrm(us, n)); 

y 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
De Moivre's Identity: 

19b (Dual number class further implementation 16g)+= (24c) <18e 

duaLnumber dual_number::pozuer(const ex & e) const { 

return dn_fromjirgjnod(argQ*e, pmv(normQ, e)); 

y 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
All algebraic routines are defined now. 

B.3. Calculation and Tests. This Subsection contains code for calculation of vari- 
ous expression. See [10] or GiNaCinfo for usage of Clifford algebra functions. 

B.3.1. Calculation of Expressions. Firstly, we output the expression of the Cayley 
transform for a generic element from subgroups N and N'. 

19c (Show expressions 19c) = (26a) 19d > 

ex XC=canonicalize_clifford((TC*X*TCI).evalm()); 

formirfaj}ut(" Cayley of the matrix x: " , XC.subs(sign=0).normalQ); 
Uses f ormula_out 19e 20e 27d. 

Then we calculate Mobius action of those matrix on a point. 

19d (Show expressions 19c) += (26a) <19c 19e> 

duaLnumber W(cliffordjnoebiusjnap(XC, P.tojnatrixQ, e).subs(sign=0).normal()), 

Wl = W.subs(lst(u=ul, v=vl)); 
formulaj)ut("Rota.tion by x: ", W); 

Uses dualjiumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b and 
f ormula_out 19e 20e 27d. 

Next we specialise the above result to the reference point. 

19e (Show expressions 19c) += (26a) <19d 19f> 

/ormu/a_OMf("Rotation of \\((u_0, v_0)\\) by \\(x\\): ", 

W.s!/bs(lst(w = uO, v = vO)).snbs(Inf = pow(y, -T)).normal().subs(y = 0).normal()); 

Defines: 

f ormula_out, used in chunks 19, 20, 23, 24a, and 28a. 

The expression for the parabolic norm. 

19f (Show expressions 19c) += (26a) < 19e 19g > 

/orw;(fo_OMf("Parabolic norm: " , P.normQ); 

Uses f ormula.out 19e 20e 27d. 

Embedding of reals into dual numbers. 

19g (Show expressions 19c) += (26a) <19f 20a > 

formulaj)ut("Real number x as a dual number: ", duaLnumber(x)); 



Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b and 
f ormula_out 19e 20e 27d. 
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The expression for the product of two points. 

20a (Show expressions 19c) += (26a) < 19g 20b > 

formnlaj}ut("Pioduct : ", P*PT); 

Uses f ormula.out 19e 20e 27d. 

The expression of the product of a point and a scalar. 

20b (Show expressions 19c) += (26a) <20a 20c > 

formulaj}ut("Piodnct by a scalar: ",a*P); 

Uses f ormula_out 19e 20e 27d. 
Expressions for the real and imaginary parts. 
20c (Show expressions 19c) += (26a) <20b 20d> 

formulaj}ut("Real part: " , P.real_partQ); 

formula j)ut(" Imag part: " , P.imag-partQ); 

Uses f ormula_out 19e 20e 27d. 

The expression for a sum of two points is too cumbersome to be printed. 

20d (Show expressions 19c) += (26a) <20c 20e> 

/ /formula_out("Add is: ", (P+Pl).normal()); 

Uses f ormula.out 19e 20e 27d. 

Linear combination of points (1, 0) and (—1,0) with coefficients a and b, for the 
linearisation presented in § 3.3.4. 
20e (Show expressions 19c) += (26a) <20d 

formalaj}ut("Lin comb of two vectors a*(l, 0)+b*(-l, 0): ", 
(a*P+b*Pl).subs(lst(u=l, v=0, ul=-\, vl=0)).normalQ); 

Defines: 

f ormula_out, used in chunks 19, 20, 23, 24a, and 28a. 

B.3.2. Checking Algebraic Identities. In this Subsection we verify basic algebraic prop- 
erties of the defined operations. 

A dual number is the sum of its real and imaginary parts. 
20f (Check identities 20f)= (26a) 20g> 

test-OUt("P is the sum Re(P) and Im(P) : ", 

P-(exJo<d\iil^umbei>(P.red^artQ)+exJo<dual^umbei>(P.imag_part()))); 

Defines: 

test.out, used in chunks 20-22, 24a, and 28a. 
Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

A dual number maid out of a real a has the norm of real part equal to a. 
20g (Check identities 20f)+= (26a) <20f 20h> 

fesLowf("The real part of a real dual number is itself: ", 
ex_to<dual_number>(dual_number(fl).reflLpflri()).Monw()-fl); 

Defines: 

test.out, used in chunks 20-22, 24a, and 28a. 
Uses dual-number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

The norm is invariant under parabolic rotations, i.e. they are in agreement with 
Defn. 3.1. 

20h (Check identities 20f)+= (26a) <20g 21a > 

fesLowf("norm is invariant: " , P.normQ-W.normQ); 



Uses test.out 20f 20g 21b 27e. 
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The product W1W2 is invariant under rotations, Prop. 3.9(iii). 

21a (Check identities 20f)+= (26a) <20h 21b > 

fesf_owf("Product is invariant: " , P*Pl.conjagateQ-W*Wl.conjugateQ}; 

Uses test.out 20f 20g 21b 27e. 

Product ww is (0, |w| 2 ),Prop. 3.9(iv). 

21b (Check identities 20f) + = (26a) <21a 21c> 

tesLowf("Product is norm squared: ", 

(P*P.conjiigateQ-dn.fromjirg_mod(ArgO, pow(P.norm(), 2)))); 

Defines: 

test .out, used in chunks 20-22, 24a, and 28a. 
Uses ArgO 16c. 

The reference point is unit under multiplication. 

21c (Check identities 20f)+= (26a) <21b 21d> 

fesf_OMf("Product \\((u, v)*(u_0, v_0)\\) is \\((u, v)\\): ",P*P0-P); 

Uses test.out 20f 20g 21b 27e. 

Addition is commutative, Prop. 3.11(i). 

21d (Check identities 20f)+= (26a) <21c 21e> 

testJ)Ut("kdd is coiranutative: ", (P+P1)-(P1+P)); 

Uses test.out 20f 20g 21b 27e. 

Addition is associative, Prop. 3.11(i). 

21e (Check identities 20f)+= (26a) <21d 21f> 

test_out(" Add is associative: ", ((P+P1)+ P2)-(P+(P1+P2))); 

Uses test.out 20f 20g 21b 27e. 

Multiplication by a scalar is commutative. 

21f (Check identities 20f}+= (26a) <21e 21g> 

test_out(" S-mult coiranutative: " , P*a-a*P); 

Uses test.out 20f 20g 21b 27e. 

Multiplication by a scalar is associative. 

21g (Check identities 20f)+= (26a) <21f 21h> 

fesLoMf("S-mult associative: " , b*P*a-a*P*b); 

Uses test.out 20f 20g 21b 27e. 

Distributive law a(wi + W2) = awi + QW2, Prop. 3.11(iv). 

21h (Check identities 20f)+= (26a) <21g 21i> 

testJ)Ut(" S-mult distributive 1: ", a*(P+Pl)-(a*P +a*Pl)); 

Uses test.out 20f 20g 21b 27e. 

Distributive law (a + b)w = aw + bw, Prop. 3.11(iv). 

21i (Check identities 20f)+= (26a) <21h 22a > 

testJ)Ut(" S-mult distributive 2: ", P*(a+b)-(P*a + P*b)); 



Uses test.out 20f 20g 21b 27e. 
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Product is commutative, Prop. 3.9(i). 

(Check identities 20f) += (26a) < 21i 22b > 

tesLowf("Product is symmetric (commutative): " , P*P1-P1*P); 

Uses test.out 20f 20g 21b 27e. 
Product is associative, Prop. 3.9(i). 

(Check identities 20f) += (26a) < 22a 22c > 

tesf.OMf("Prod is associative: ", (P*P1)*P2-P*(P1*P2)); 

Uses test.out 20f 20g 21b 27e. 

Product and addition are distributive, Prop. 3.11(ii). 

(Check identities 20f}+= (26a) <22b 

fesLoMf("Product is distributive: " , (P+P1)*P2-(P*P2+P1*P2)); 



Uses test.out 20f 20g 21b 27e. 

B.4. Induced Representations. Here we calculate the basic formulae for Section 4. 

B.4.1. Encoded formulae. This routine encodes the map s : M 2 — > SL 2 (R) (4.1). 

(Induced representations routines 22d) = (24e) 22e > 

ex s_map(const ex & u, const ex & v) { 
return matrix(2, 2, lst(u,u,0,l)); 

} 

ex s_nwp(const ex & P) { 
if (P.nopsQ = 2) 

return sjnap(P.op(0), P.op(l)); 
cerr "s_map() error: parameter should have two operands" <C endl; 
return s_map(P,l); 

} 



This routine encodes the map r : SL 2 (M) — > H (4.2). The first parameter is an 
element of SL 2 (R), the second — is a generic element of subgroup H. 

(Induced representations routines 22d) += (24e) < 22d 22f > 

ex r_map(const ex & M, const ex & K) { 
ex Kl=K.evalmQ, K2; 

1st vars = (z's_fl<symbol>(Ja.op(2)) ? lst(X2.op(2)) : lst(Kl.op(l))); 
if (/s_a<symbol>(KLop(3))) { 
vars = vars.append(Kl.op(3)); 

K2 = Kl.subs(lsolveQst((M*Kl).evalmQ.op(2)=0), vars)).subs(Kl.op(3)=l); 
} else 

K2 = Kl.subs(lsolve(lst((M*Kl).evalm().op(2)=0), vars)); 
return pow(K2, -l).evalm(); 

} 



This is the inverse s 1 of the above map s. 

(Induced representations routines 22d) += (24e) < 22e 23a > 

ex sJnvfconst ex & M, const ex & K) { 

ex MK=(M*pow(r_map(M,K),-l)).evalm(); 
ex D=MK.op(3).subs(x=l).normal(); 

return matrix(l, 2, lst((MK.op(l).subs(x=l).normal()^D).normal(), 

(MK.op(0).subs(x=l).normal()^D).normal())); 

} 
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This is a matrix form of the above inverse map sJnvQ. 

23a (Induced representations routines 22d) += (24e) < 22f 

ex sjnvjn{const ex & M, const ex & K) { 
return (M*pow(r_map(M,K),-l)).evalmQ; 

} 



B.4.2. Caculation of induced representation formulae. Firstly we define a generic ele- 
ment M of SL 2 (R). 

23b (Induced representations 23b) = (26e) 23c> 

ex M=matrix(2,2, \st{a,b,c,d)), H; 

We consider the three cases. 

23c (Induced representations 23b) += (26e) <23b 23d > 

char* cflses[]={"Elliptic", "Parabolic (\\ (N~\\prime\\) ) ", "Hyperbolic"}; 

Defines: 

cases, used in chunk 23e. 

In the those cases subgroups holds a generic element of a subgroup H, see (2.2) 
and (2.4). 

23d (Induced representations 23b) += (26e) <23c 23e> 

ex swbs>roMps=lst(matrix(2, 2, \st(x,-y,y,x)), 
matrix(2, 2, lst(l,0,y,l)), 
matrix(2, 2, lst(x,y,y,x))); 

Now we run a cycle over the three cases. . . 
23e (Induced representations 23b) += (26e) <23d 23f> 

for(int i=0; f<3; ;'++){ 
H=subgroups[i]; 

cout <C cases[i] <C " case of induced representations\\\\" <C endl; 
//formula_out("M*H: ", (M*H).evalm()); 

Uses cases 23c and f ormula_out 19e 20e 27d. 

. . . and output expression of r (4.2),. . . 

23f (Induced representations 23b) += (26e) <23e 23g> 

formula .out {"map \\(r(M)\\): " , r_map(M,H)); 

Uses f ormula_out 19e 20e 27d. 

. . . matrix form of the inverse s _1 (4.2),. . . 

23g (Induced representations 23b) += (26e) <23f 23h> 

formula j)ut(" map \\(s"{-l>(M)\\) : ", sJnv_m(M,H).subs(a=(l+b*c)~d).normal()); 

Uses f ormula.out 19e 20e 27d. 

. . . expression for the argument of the character in (4.4),. . . 

23h (Induced representations 23b) += (26e) <23g 24a > 

formulajout(" character : ", rjnap(M*sjnap(P),H)); 



Uses f ormula.out 19e 20e 27d. 



24 



VLADIMIR V.KISIL 



. . . and finally the action (4.3) of SL 2 (R) on the homogeneous space. 

24a (Induced representations 23b) += (26e) <23h 

/ormwZfl_owf("Moebius map: ", sJnv(M*sjnap(P.tojnatrix()),H)); 

fesf_OMi("Moebius map is given by the imaginary unit: ",sJnv(M*sjnap(P),H)' 

Clifford jnoebiusjnap{a*one, b*eO, -c*eO, d*one,P.tojnatrix(),e).subs(sign=i-l)); 
cout <C (latexout ? "\\vspace{2mm}\\hrule" : 

ii ii ) ^ m( n- 

y 

Uses f ormula.out 19e 20e 27d and test.out 20f 20g 21b 27e. 

B.5. Program Outline. Here is the outline how we use the above parts. 
Routines for dual_number are collected in a separate library. We start from the 
definition dual_number class in the header file. 

24b (duamum.h24b) = 
(Initialisation 24f) 
(N-Nprime separation 16a) 
(Dual number class declaration 28b) 
(Additional routines declarations 28a) 

Here is the file with the implementation. 

24c (dualnum.cpp 24c) = 

#incfude <dualnum.h> 
(Algebraic procedures 16d) 
(Dual number class implementation 29c) 
(Dual number class further implementation 16g) 
(Output routines 27d) 

B.5.1. Test program outline. Firstly we load dual_number support. 

24d (*24d)= 24e> 
#include <dualnum.h> 

The rest of the program makes all checks. 

24e (*24d)+= <24d 
(Definition of variables 25a) 
(Test routine 25f) 

(Induced representations routines 22d) 
(Main procedure 26b) 

This is the initialisation part 

24f (Initialisation 24f)= (24b) 
#include <cycle.h> 
#include <fstream> 
using namespace std; 
using namespace GiNaC; 



#include "ginac-utils .h" 
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B.5.2. Variables. These realsymbols are used in our calculations. 
25a (Definition of variables 25a) = (24e) 25b > 

realsymbol u("u"), v("v"), ul("u"'), vl("v"'), m2("u' ' "), v2("v' ' "), 
fl("a"), fo("b"), c("c"), d("d"), x("x"), y("y"), 

Finally this variable keeps the signature of the metric space. 

25b (Definition of variables 25a) += (24e) < 25a 25c > 

sign("s", "Wsigma"); 

This an index used for the definition of Clifford units. 

25c (Definition of variables 25a) += (24e) <25b 25d> 

varidx m«(symbol("mu", "\\mu"), 2); 

Three generic points which are used in calculations. 

25d (Definition of variables 25a) += (24e) <25c 25e> 

dual-number P(u, v), Pl(ul, vl), P2(u2, v2); 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
Here we define a parabolic Clifford units eO, el. 

25e (Definition of variables 25a) += (24e) <25d 

ex e = Clifford Mnit(mu, diagjnatrix(lst(-l, sign))), 
eO = e.subs(mn=0), 
el = e.subs(mn=l), 
one = diracJDNEQ; 



B.5.3. Test routine. This routine make the same sequence of checks for both cases 

of subgroups N and N'. 

First we define the reference point (uO,vO). 

25f (Test routine 25f)= (24e) 25g> 

void parab_rot^ub() { 
exX, 

uO=ArgO, vO=v.fromjiorm{uO, 1), 
PO=matrix(l, 2, lst( uO, vO)), 



Defines: 

parab_rot_sub, used in chunk 26d. 
Uses ArgO 16c. 

These two matrices define the Cayley transform and its inverse. 

25g (Test routine 25f)+= (24e) <25f 25hi> 

TC=matrix(2, 2, lst(owe, -el, -el, one)), 
TCI=matrix(2, 2, lst(one, el, el, one)); 



For the subgroup N we consider upper- triangular matrices, for N' — lower- triangular. 

25h (Test routine 25f)+= (24e) <25g 26a > 

if (is sub group _N) 

X=matrix(2, 2, lst(one, eO*x, 0, one)); 

else 

X=matrix(2, 2, lst(o«e, 0, e0*x, one)); 
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Common part of test routine. 

(Test routine 25f)+= (24e) <25h 

(Show expressions 19c) 
(Check identities 20f) 

cout <C (latexout ? "\\vspace{2mm}\\hrule" : 

" " ) < endl; 

} 



B.5.4. Main procedure. It just calls the test routine, calculates the induced represen- 
tation and draws a few pictures. 
We output formulae in KTrjX mode. 

(Main procedure 26b) = (24e) 26c> 

int main(){ 

latexont=true; 



Defines: 

main, never used. 

Preparation of output stream. 

(Main procedure 26b) += (24e) <26b 26d> 

cout <S boolalpha; 
if (latexout) 

cout <C latex; 



Now we call the test routine for both N and N' subgroups. 

(Main procedure 26b) += (24e) <26c 26e> 

issubgroupJSl = true; 
parabsotsubQ; 
issubgroup-N = false; 
parabsotsubQ; 

Uses parab_rot_sub 25f. 

Calculation of induced representations formulae. 

(Main procedure 26b) += (24e) <26d 26f> 

(Induced representations 23b) 

And we finishing by drawing several pictures for Figs. 1 and 2. 

(Main procedure 26b) += (24e) <26e 

(Drawing pictures 26g) 
} 



B.6. Drawing Orbits. To draw cycles we use cycle library [12]. 
Elliptic orbits (circles). 

(Drawing pictures 26g)= (26f) 27a > 

ofstream asymptoie("parab-rot-data. asy"); 
asymptote <C "path [] K="; 
for(int i=0; i<6; i++) 

cycle2D(lst(0,0),e.sufos(s!,gw=-l),i*i*.04) 

.asyjpath(asymptote, -1.5, 1.5, -2, 2, 0, (z>0)); 
asymptote endl; 
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Hyperbolic orbits. 

27a (Drawing pictures 26g)+= (26f) <26g 27b > 

asymptote <C "path [] A="; 
for(int i=0; i<6; 

cycle2D(lst(0,0),e.sMfe(s!gn=l),-i«*.04) 

.asyjpath(asymptote, -1.5, 1.5, -1.5, 2, 0, 0>0)); 
asymptote <c " ; " <c end/; 



Parabolic orbits, subgroup N . 

27b (Drawing pictures 26g)+= (26f) <27a 27c> 

asymptote <C "path [] N="; 
for(int ;=0; /<6; 

cycle2D(l,lst(0,numeric(l,2)),numeric(i,2)-l,e.sMfos(s/^«=0)) 
.asyjpatli(asymptote, -1.5, 1.5, -2, 2, 0, (/>0)); 
asymptote <C " ; " <C end/; 



Parabolic orbits, subgroup N'. 

27c (Drawing pictures 26g)+= (26f) <27b 

asymptote <C "path[] Nl="; 
for(int i=0; i<5; i++) 

cycle2D(.5*;*/*f+l,lst(0,numeric(l,2)),-l,e.s!(fos(sf^n=0)) 
.asyjpath(asymptote, -1.5, 1.5, -1.5, 2, 0, (z>0)); 
asymptote <C " ; " <C end/; 



asymptote. closeQ; 



B.6.1. Output routines. We use standardised routines to output results of calcula- 
tions. 

27d (Output routines 27d) = (24c) 27e> 

void formula _ou/(char* S, ex F) { 

com/ < S <C (latexout ? "\\ (":"")« F <C (latexout ? "\\)\\\\" : "" ) 
< end/; 

} 



Defines: 

f ormula_out, used in chunks 19, 20, 23, 24a, and 28a. 

This routine is used to check identities. 

27e (Output routines 27d)+= (24c) <27d 

void test _owf(char* S, ex T) { 

cout < S <C (latexout ? "Wtextbf {" : "*" ) 

<C (z's_fl<dual_number>(T) ? ex_fo<dual_number>(T).Mormfl/()j's_zero() : 

T.evalmQ.normalQ.is-zerojnatrixQ) <C (latexout ? ">\\\\" : "*" ) 
< end/; 

} 



Defines: 

test.out, used in chunks 20-22, 24a, and 28a. 
Uses dual.number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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Here is declarations of additional routines for the header file. 

28a (Additional routines declarations 28a) = (24b) 
bool latexout; 

duaLnumber dn_fromjirgjnod(const ex & a, const ex & n); 
ex o_/rom_norm(const ex & it, const ex & n); 
void test xut (char* S, ex T); 
void/on«wZfl_o«i(char* S, ex F); 

Uses dual .number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b, 
f ormula_out 19e 20e 27d, and test.out 20f 20g 21b 27e. 

B.7. Header and Implementation of the dual .number Class. 

B.7.1. Header File for dualjiumber. We use the standard GiNaC machinery do de- 
fine dual .numbers as derived of the class basic. 
28b (Dual number class declaration 28b) = (24b) 28c > 

class duaLnumber : public basic 

{ 

GINAC.DECLAREJlEGISTERED_CLASS(dvLal-riumber, basic) 
static const tinfostaticj return JypeJinfo_static[256]; 

Uses dual .number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b and 
tlnf o.static.t 29c. 

Public methods (constructors, algebraic, etc.) 

28c (Dual number class declaration 28b) += (24b) <28b 28d> 

public: 

(Public methods 14) 
(Technical methods 28f) 

We redefine protected methods for printing only. 

28d (Dual number class declaration 28b) += (24b) <28c 28e> 

protected: 

void dojprin t(const print xontext & c, unsigned level) const; 
void dojprintJatex(const print latex & c, unsigned level) const; 

Private members: two components of a duaLnumber. 

28e (Dual number class declaration 28b) += (24b) <28d 29a > 

protected: 

ex uxomp; 
ex vxomp; 

>; 

The following methods are needed for GiNaC to work properly. 

28f (Technical methods 28f)= (28c) 

dualjiumber normalQ const { return dual_number(i;_comp. normalQ, vxomp. normalQ); } 

dualjiumber si;fes(const ex & e, unsigned options = 0) const; 

bool is-zeroQ const; 

bool isxqual(const ex & other) const; 

size J nopsQ const { return 2; } 

ex op(sizeJ i) const; 

ex & letxp(sizeJ i); 



Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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We overload standard algebraic operations for dual_number. 

29a (Dual number class declaration 28b) += (24b) <i28e 29b > 

const dual_number operator+(const dual_number & Ih, const duaLnumber & rh); 
const dual_number operator-(const dual_number & Ih, const dual_number & rh); 
const dual_number operator*(const dual_number & Ih, const duaLnumber & rh); 
const dual_number operator*(const dual_number & Ih, const ex & rh); 
const dual_number operator*(const ex & Ih, const dual_number & rh); 
const dual_number operator^(const dual_number & Ih, const ex & rh); 
const dual_number operator^(const ex & Ih, const dual_number & rh); 
const dual_number operator^(const dual_number & Ih, const ex & rh); 



Defines: 

dual-number, used in chunks 14—20, 25d, 27-30, and 33. 

29b (Dual number class declaration 28b) += (24b) < 29a 

dual-number dn_fromjirgjnod(const ex & a, const ex & n); 
1 1 End of "header" 



Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

B.7.2. Standard Implementation Part. The implementation uses standard GiNaC tech- 
nique. 

29c (Dual number class implementation 29c) = (24c) 29d > 

GINACJMPLEMENTJLEGISTEREDJ0LASS-OPT(d\i*l_n\imber, basic, 

print _func<printxontext>(lkdud\jcmmber::dojprint). 
print _func<printJatex>(&dud\jcmmber::dojprint latex)) 

const HnfastaticJ dual_number::refwniJype_h'n/o_staf!c[256] = {{}}; 

DEFA ULT .ARCHIVING (dual_number) 

Defines: 

dual-number, used in chunks 14—20, 25d, 27-30, and 33. 
tinf o_static_t, used in chunk 28b. 

B.7.3. Implementation of Constructors. Default constructor. 

29d (Dual number class implementation 29c) += (24c) <29c 29e> 

duaLnumber: :dual_number() : inherited(&cdual^mmbet::tinfo^tatic), uxomp(Q), vjcomp(0) 
{ 

setflag(status.fl.ags::not shareable); 

y 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

Constructor from two components. 

29e (Dual number class implementation 29c) += (24c) < 29d 30a > 

duaLnumber: :dual_number(const ex & a, const ex & b) : inherited(&cd\ial-number::tinfostatic), 

uxomp(a), vjcomp(b) 

{ 

} 



Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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Constructor from a single expression. It may contain two components. . . 

30a (Dual number class implementation 29c) += (24c) <29e 30b > 

dual_number::dual_number(const ex & P) : inherited(&d\ial-t\umber::tinfo-Static) 
{ 

if (isji<lst>(P) V z's_fl<matrix>(P) V 7s_a<dual-number>(P)) { 
uxomp = P.op(0); 
vxomp = P.op(l); 

Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

... if it is a real expression we embed it into duaLnumber. . . 

30b (Dual number class implementation 29c) += (24c) < 30a 30c > 

} else if (P " .imag_partQ.isjzero{)) { 
uxonvp = ArgO; 

vxomp = v_fromjiorm{ArgO, P); 
Uses ArgO 16c. 

... or if its a complex expression we decompose it into the real and imaginary parts. 

30c (Dual number class implementation 29c) += (24c) <30b 30d> 

> else { 

uxomp = P.realjpartQ; 
vxomp = P .imagjpartQ; 

} 

} 

Comparison routine. 

30d (Dual number class implementation 29c) += (24c) <30c 30e> 

int dual_number: -.compare same Jype{const basic & other) const 
{ 

GINACJLSSERT(isji<dual-B,wnbei>(otherj); 

const duaLnumber &o = static_cast<const dual_number &>(other); 

int cmpval = uxomp. compare(o.op(Q)); 
if (cmpval^O) return cmpval; 

return vxomp. compare(o.op(l)); 



Defines: 

dualjiumber, used in chunks 14-20, 25d, 27-30, and 33. 
Equality of two dual numbers. 
30e (Dual number class implementation 29c) += (24c) <30d 31a > 

bool dual_number::!s_eqwfl/(const ex & other) const 
{ 

GINACJLSSERT(isM<dualJxambei>(other)); 

const duaLnumber &o = static_cast<const dual_number &c>(other); 
return uxomp.isxqual(o.op(0)) A vxomp. is xqual(o.op(l)); 



Uses duaLnumber 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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B.7.4. Printing. Default printing. 

31a (Dual number class implementation 29c) += (24c) <30e 31b > 

void dual_number::do_pnnf(const printxontext & c, unsigned level) const 
{ 

c.s < " ["; 

ujcomp.print(c); 

c.s< 

vxomp.print(c); 
c.s < "] "; 

} 

Defines: 

dualjmmber, used in chunks 14—20, 25d, 27-30, and 33. 

KTthX printing. 

31b (Dual number class implementation 29c) += (24c) < 31a 31c > 

void duaLnumber::do_praif Jata(const print Jatex & c, unsigned level) const 
{ 

c.s < "\\left(\\begin{array}{cc}"; 
ujcomp.print(c); 

c.s <C 

w_comp.prmf(c); 

c.s <C "\\end{array}\\right) "; 

} 

Defines: 

dualjiumber, used in chunks 14—20, 25d, 27-30, and 33. 

B.7.5. Overloading algebraic operations. Addition. 

31c (Dual number class implementation 29c) += (24c) < 31b 32a > 

const dual_number operator+(const duaLnumber & Ih, const duaLnumber & rh) 
{ 

return lh.add(rh); 

y 

const dual_number operator-(const dual_number & Ih, const dual_number & rh) 

{ 

return lh.sub(rh); 

y 



Defines: 

dualjiumber, used in chunks 14-20, 25d, 27-30, and 33. 
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Multiplication. 

32a (Dual number class implementation 29c) += (24c) < 31c 32b > 

const dual_number operator*(const dual_number & Ih, const ex & rh) 
{ 

return lh.mul(rh); 

} 

const dual_number operator*(const ex & Ih, const dual_number & rh) 

{ 

return rh.mulQh); 

y 

const dual_number operator*(const dual_number & Ih, const duaLnumber & rh) 

{ 

return lh.mnl(rh); 

y 

Defines: 

dualjiumber, used in chunks 14-20, 25d, 27-30, and 33. 
Division. 

32b (Dual number class implementation 29c) += (24c) < 32a 33a > 

const dual_number operators (const dual_number & Ih, const dual_number & rh) 

{ 

return lh.mul(rh.power(-l)); 

y 

const dual_number operators (const ex & Ih, const dual_number & rh) 

{ 

return rh.power(-l)*lh; 

y 

const dual_number operators (const dual_number & Ih, const ex & rh) 
{ 

return lh.mul(pow(rh, -1)); 

y 



Defines: 

dualjiumber, used in chunks 14-20, 25d, 27-30, and 33. 
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B.7.6. Component-related functions. 

33a (Dual number class implementation 29c) += (24c) < 32b 33b > 

ex dual_number: :op(size J i) const 
{ 

GINAC-ASSERT(i<nopsQ); 

switch (i) { 
case 0: 

return iijcomp; 
case 1: 

return vxomp; 
default: 

thiow(std:\invalidjirgument("d\ial_nxmber : :op(): requested" 

" operand out of the range (2)")); 

} 

} 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

33b (Dual number class implementation 29c) += (24c) < 33a 33c > 

dual_number dual_number::sufos(const ex & e, unsigned options) const 
{ 

exmap m; 

if (e.info(info.flags::list)) { 
1st / = ex.fo<lst>(e); 

for (lstv.const Iterator i = l.beginQ; i ^ I.endQ; ++i) 
m.insert(std::make_pair(i->op(0), i— >op(l))); 
} else if (;s_a<relational>(e)) { 

m.insert(std::makejpair{e.op(Q), e.op(\))); 
} else 

thiow(std:jnvalidjirgiiment("dual_n.\mber : :subs() : the parameter" 

" should be a relational or a 1st")); 

return exJo<d\ial-number>(inherited::subs(m, options)); 
} 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 

33c (Dual number class implementation 29c) += (24c) < 33b 

ex & dual_number::/ef_op(size_f i) 
{ 

GINAC.ASSERT(i<nopsQ); 

ensured/ jnodifiableQ; 
switch (i) { 
case 0: 

return ujzomp; 
case 1: 

return vxomp; 
default: 

throw(std:\invalidjjrgument("dxial_n.\mber : :let_op() : requested operand" 

" out of the range (2)")); 

} 

} 

Uses dual_number 29a 29a 29a 29a 29a 29a 29a 29a 29c 30d 31a 31b 31c 31c 32a 32a 32a 32b 32b 32b. 
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This is the end of dual_number implementation. 
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